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ABSTRACT 

Approximately 1% of low redshift (z < 0.3) optically-selected type 2 Active 
Galactic Nuclei (AGNs) show a double-peaked [Om] narrow emission line profile 
in their spatially-integrated spectra. Such features are usually interpreted as due 
either to kinematics, such as biconical outflows and/or disk rotation of the narrow 
line region (NLR) around single black holes, or to the relative motion of two distinct 
NLRs in a merging pair of AGNs. Here we report follow-up near infrared (NIR) 
imaging and optical slit spectroscopy of 31 double-peaked [Om] type 2 AG Ns drawn 



from th e Sloan Digital Sky Survey (SDSS) parent sample presented in iLiu et al. 



(|2010bh . The NIR imaging traces the old stellar population in each galaxy, while the 
optical slit spectroscopy traces the NLR gas. These data reveal a mixture of origins 
for the double-peaked feature. Roughly 10% of our objects are best explained by 
binary AGNs at (projected) kpc-scale separations, where two stellar components 
with spatially coincident NLRs are seen. ~ 50% of our objects have [Om] emission 
offset by a few kpc, corresponding to the two velocity components seen in the SDSS 
spectra, but there are no corresponding double stellar components seen in the NIR 
imaging. For those objects with sufficiently high quality slit spectra, we see velocity 
and/or velocity dispersion gradients in [Om] emission, suggestive of the kinematic 
signatures of a single NLR. The remaining ~ 40% of our objects are ambiguous, 
and will need higher spatial resolution observations to distinguish between the two 
scenarios. Our observations therefore favor the kinematics scenario with a single 
AGN for the majority of these double-peaked [Om] type 2 AGNs. We emphasize 
the importance of combining imaging and slit spectroscopy in identifying kpc-scale 
binary AGNs, i.e., in no cases does one of these alone allow an unambiguous 
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identification. We estimate that ~ 0.5-2.5% of the z < 0.3 type 2 AGNs are kpc-scale 
binary AGNs of comparable luminosities, with a relative orbital velocity > 150kms"'. 

Based in part on observations obtained with the 6.5 m Magellan telescopes located 
at Las Campanas Observatory, Chile, and with the Apache Point Observatory 3.5 m 
telescope, which is owned and operated by the Astrophysical Research Consortium. 

Subject headings: black hole physics — galaxies: active — galaxies: interactions — 
galaxies: nuclei — galaxies: Seyfert — quasars: general 



Introduction 



Binary supermassive black holes (SMBHs) have long been proposed theoretically (e.g. jBegelman et al 



1980h . Galaxy mergers frequently occur and it is now widely appreciated that most bulge-dominated 
merging galaxies harbor central SMBHs. Thus the formation of binary SMBHs from galaxy merg- 
ers seems inevitable. Major mergers between two gas-rich galaxies of comparable mass also pro- 
vide a viable means to c hannel a large am ount of gas into the central region to trigger starbursts 
and quasar activity (e.g., lHernquistl[l989 ^). On the other hand, observationally confirmed binary 
SMBHs are surprisingly scarce. Abo ut 0.1% of quasars are in physical pairs with p rojected sep- 
arations of tens to hundreds of kpc (IHennawi et al.l l2006l 1201 Ol : iMyers et al.l l2008h . and only a 
handful o f confirmed binaries with projected separ ations on kpc scales are known in the litera- 
ture (e.g., iKomossa et al.ll2003l : iBianchi et al.ll20() 8 [). On pc scales, the re is only one confirmed 
binary SMBH separated by ~ 7pc in projection ( Rodriguez et al. 20061) . detected with the Very 
Long Baseline Array (VLBA). There are essentially n o confirmed sub-pc b inaries, although several 
candidates have been pro posed, such as OJ287 (e.g.. IValtonen et al.ll2008h and SDSSJ1536-I-0441 
(IBoroson & Laueiil2009|) . for both of which alternative, non-binary interpretations exist. The fre- 
quency of observed SMBH binaries is thus far below the expectations from the observed (and 
simulation-predicted) galaxy merger rate, and is in tension with the hypothesis that quasar activity 
is triggered by major mergers. 

If mergers do produce binary SMBHs, then the most natural explanation for the apparent 
deficit of SMBH binaries is that only a tiny fraction of them are detectable; both BHs must be 
active for a detection and either their spatial separations must be resolvable, or orbital motion 
can be unambiguously inferred. So far the searches for binary AGNs/quasars with separations 
below several kpc are highly incomplete due to the stringent spatial resolution requirement. If 
quasar activity only occurs at the late stage of a merger (as simulations suggest), then most binaries 
will be difficult to detect. A further complication is that at relatively low AGN luminosity (i.e., 
^boi ^ 10"^^ erg s"^), mergers may not be necessary to trigger BH activity, and alternative fueling 
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routes (secular processes) may suffice. Simulations predict that the effects of mergers on BH 
accre tion start to become important at late stages when a pair of BHs is separated by a few kpc 
(e.g. jHopkins et al.ll2005[) : this is a scale of which one can still resolve the two nuclei with current 
instruments. It is therefore important to quantify the fraction of kpc-scale binary AGNs. This 
fraction will shed light on the AGN triggering mechanism, and individual binaries will provide 
ideal test beds to understand the interplay between AGNs and hosts in these merging systems. 

The most efficient and systematic way to find these kpc-scale binary AGNs is an imaging 
survey conducted in the optical/NIR with sub-arcsec resolution (to find kpc separati on pairs) and 
the ability to identify AGNs. The seeing in imaging data from the SDSS flYork et al.ll200Q) is only 
~ 1.4", but SDSS does recover some very low-redshift kpc-scale binary AGNs (Liu et al. 2010c, 
in preparation), where the two AGNs are separated by more than a few arcsec and were observed 
spectroscopically with separate fibers. At higher redshift, such kpc-scale binary AGNs will fall 
within the 3" fiber and spatial information within the fiber is completely lost. A different approach 
is to select candidates with spectrosco pic features t hat m ay hi nt at a binary and follow them up 



Heckman et al.lll98ll : IZhou et al. 



to confirm their binary AGN nature. iGerke et al.l (|2007[) and iComerford et al.l (|2009aD revived 
the idea of selecting candidate kpc-sca le binary AGNj ^ based on double-peaked [Om] lines (e.g.. 



2004|) . where the hope is that the two velocity components of the 



[Oiii] line originate from two distinct narrow line regions (NLRs), each associated with its own 
BH. These authors found that in two objects at z ~ 0.7 the two [Om] components were spatially 
offset by a few kpc in thei r slit spectra, and claime d them to be kpc-scale binary AGNs. A third 
candidate was reported by IComerford et al.l (l2009b|) based on spatially of fset [Onil ernission in a 
slit spectrum and a pair of resolved double nuclei in HST imaging (but see lCivano et al.ll2010l for 
an alternative interperation for this system). 

Several independent searches for AGNs w ith double-peaked [O m] lines have sinc e been con 



ducted using the SDSS spectroscopic database: iLiu et al.1 (1201 Obi hereafter Paper I) and lWang et al. 



(120091) selected such objects from type 2 AGNs, while ISmith et al.l (|2010l) used a sample of type 1 
AGNs (mixed with some type 2 objects due to SDSS pipeline misclassification). More than 200 
double-peaked [Om] AGNs were identified (based on visual inspection of SDSS spectra) with a 
median redshift z ~ 0. 15, amounting to about 1% of the total AGN population. With very few 
excepti ons, these objects do not appear to have double nuclei in their optical SDSS images (e.g.. 
Paper I: ISmith et al.ll2010|) : and of course the SDSS fiber spectra do not indicate if the two velocity 
components are spatially offset by a few kpc. The bulk properties of the double-peaked AGNs are 
not dramatically different from those of the normal AGN population (Paper I). 



Not all of these double-peaked [Om] objects are kpc-scale binary AGNs. Spatially resolved 
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studies of nearby Seyferts frequently show a quite dynamic picture of the NLR, involving outflows, 
inflows, rotation, and interactions with a radio jet. Indeed, double-peaked [O nil features were 
already seen in early NLR studies (e.g.. lSargentill972l : lHeckman et al.lll981uVeilleuxlll991|) . One 
famous example is Mrk 78, which has a bipolar outflowing NLR structure, broadly ali gned with the 
linear radio jet (e.g. Jwhittle & Wilsonll2004l : Iwhittle et alJllOOsI: Ipischer et alJboiol): its soatiallv 



integ rated spectrum shows a clear double-peaked [Om] line profile (see fig. 1 of iHeckman et al. 



198lk Spatially resolved studies with HSThave revealed more such cases (such as NGC 1068 and 



NGC 4151, e.g. JCrenshaw & Kraemedl200(]l : ICrenshaw et al.ll200(j : IVeilleux et albOOlUPas et al. 



20051) . where the NLR is dominated by outflows and/or rotation. For these dynamic NLRs, double- 
peaked features will naturally arise in spatially-integrated spectra when viewed at proper angles. 
Moreover, these outflows or rotating gaseous disks can extend as far as a few kpc. Thus detecting 
two [Oiii] peaks offset by a few kpc in slit spectra alone is not sufficient to conclude that the object 
is a kpc-scale binary AGN. All these complications simply reflect the diverse nature of NLR gas 
dynamics, and follow-up observations are needed to identify bona fide kpc-scale binary AGNs 
from the double-peaked [O m] sample. 

To that end, we have been conducting follow-up observations of our double-peaked [Om] 
sample (Paper I) with ground-based NIR imaging and slit spectroscopy. The NIR imaging offers 
spatial resolution ~ 0.6" (much better than the typical ~ 1.4" seeing for optical SDSS imaging), 
and probes the old stellar populations (mostly stars in the bulge). Unlike the NLR gas, the stellar 
populations are unlikely to be involved in outflows, thus are a much less ambiguous indicator of the 
binary nature than using the NLR gas as a tracer. Slit spectroscopy probes the spatial distribution 
of NLR gas emission, giving the spatial information unavailable in the SDSS fiber spectra. Once 
a pair of nuclei are resolved in the NIR imaging, slit spectroscopy is required to associate the two 
[Om] peaks with each of the two nuclei. If, as is the case for Mrk 78, the NIR image shows 
a smooth stellar distribution while the optical spectrum shows two spatially offset [Om] peaks 
without NIR counterparts, then the double-peaked feature seen in the spatially integrated spectrum 
is more likely due to the kinematics of the NLR of a single AGN. Our imaging and spectroscop y 
follow-up observations have led to the discovery of several candidate binaries (|Liu et al.ll2010a|) . 
In these objects, we see spatially resolved double nuclei in the NIR images, whose locations are 
coincident with the two [Om] components in the slit spectra. The two [Om] components have 
a relative velocity offset of a few hundred kms"' , which leads to the double-peaked profile in the 
SDSS spectrum. The two component s in each case have em ission line flux ratios that place them in 
the AGN region of the BPT diagram ([Baldwin et al.lll98lb . These cases strongly suggest that they 
are binary AGNs at ~ kpc (projected) separations, which are co-rotating along with their stellar 
bulges and are ionizing their individual NLRs in a merging pair of galaxies. 



Recently, the interest of these double-peaked [Oii i] objects has led to two NIR imagin g stud- 
ies with ground-based adaptive optics (AO) systems (|Fu et al.ll2010l : iRosario et al.ll201 1|) . With 
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diffraction-limited resolution (~ 0.1"), these studies were able to resolve even closer double nu- 
cleus in the NIR than our imaging data under natural seeing. However, optical spectroscopy is 
still needed to register the resolved double nucleus with the [Oiii] emission in order to confirm the 
binary AGN nature, which becomes challenging on these smallest scales from the ground. 

In this paper we present a summary of our imaging and spectroscopy follow-up for 3 1 objects 
in our parent sample of double-peaked [Oiii] type 2 AGNs presented in Paper I. The details of the 
parent sample (167 objects) can be found in that paper. Here we briefly describe the construction 
of this sample. We started from a sample of ~ 15,000 type 2 AGNs mostly from the MPA- 
JHU SDSS DRV galaxy sample^ with the following criteria: (1) the rest-frame wavelength ranges 
[4700, 5100]A and [4982, 5035] A centered on the [Oiii] A5007 line have median signal-to-noise 
ratio (S/N) > 5 pixel"^ and bad pixel fraction < 30%; (2) the [Oiii] A5007 line is detected at > 5a 
and has a rest-frame equivalent width (EW) > 4A; 3) the line flux ratio [Oiii] A5007/H/3> 3 if 
z > 0.33, or the diagnostic line ratios [Oiii] A5007/H/3and [Nii] A6584/Hq; lie above the theoretical 
upper limits for star-formation excitation from Kewley et al. (2001) on the BPT diagram (Baldwin 
et al. 1981) if z < 0.33. All AGNs were then visually inspected and those with well-detected double 
peaks in both [Oiii] A4959 and [Oiii] A5007 with similar profiles were included in our final sample. 
We did not include those with complex line profiles such as lumpy, winged, or multi-component 
features. 

The structure of the paper is as follows. We describe our observations and data reduction in 
^ followed by discussions of individual objects in ^ where we categorize objects into binaries 
( ^3.11) , NLR kinematics around single AGNs ( S3. 21 ). and ambiguous cases ( S3. 31) . We discuss the 
frequency of kpc-scale binary AGNs in ® and conclude in ^ Throughout this paper we adopt a 
flat ACDM cosmology with Qq = 0.3, = 0.7 and Hq = 70kms-^Mpc-^ 



2. Observations and Data Reduction 



2.1. NIR Imaging 



We obtained ^^-band (or 7-band if is unavailable) images for 61 objects in our double- 
peaked sample during si x nights in two obs erving runs using the Persson's Auxiliary Nasmyth 
Infrared Camera (PANIC; Martini et al.ll2004 ) on the 6.5 m Magellan I (Baade) telescope. Here we 
reporj^ 31 of these objects for which we also have slit spectroscopy, and will publish the remaining 
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imaging data when slit spectroscopy data become available. The first observing run was on the 
nights of 2009 December 29 through 2010 January 2 UT and the second was on the night of 
2010 May 30 UT. The typical observing procedure consisted of two sequences each dithered at 
nine positions with 25" offsets. We obtained four 15-second exposures at each p osition so that 
the ty pical total exposure time per target was 18 min. We observed standard stars (I Persson et al. 



19981) at the beginning, middle, and end of each night. The observing conditions were clear but not 



photometric, with seeing ranging between 0".5 and 0".8 in the optical (through an RG610 filter) 
during the first run and 0".7 and 1".5 during the second. Table [T] lists total exposure times and 
seeing measured from field stars whenever available or from adjacent observations otherwise. 

PANIC has a 2' X 2' field of view (FOV) and 0" 125 p ixels. We reduced PANIC data using the 
Carnegie Supernova Project pipeline ( Hamuy et al. 20061) following standard procedures, includ- 
ing spatial-distortion correction, dark subtraction, bad-pixel masking, flat fielding (using twilight 
flats), sky subtraction, and aligning and st acking of the dithered frames. We determined and 
J photometric zero points using 2MASS dSkrutskie et al.l 120061) magnitudes of field stars when 
available, or the standard stars we observed during each night. 



2.2. Optical Slit Spectroscopy 

We conducted optical slit spectroscopy using the Low-Dispersion Survey Spectrograph (LDSS3) 
on the 6.5 m Magellan II (Clay) telescope and the Dual Imaging Spectrograph (DIS) on the Apache 
Point Observatory 3.5 m telescope. To date we have observed 31 objects with NIR imaging data. 
All objects which appear to have double nuclei in their NIR images were observed with slit spec- 
troscopy. 

Our LDSS3 runs were on the nights of 2010 January 12 through 14 UT. The observing con- 
ditions were clear but not photometric, with seeing ranging between 0".7 and I'M. LDSS3 has a 
8'. 3 diameter FOV and 0".188 pixels. We employed a l"x4' long-slit with the VPH-Blue grism 
to cover H/3 and [Om] in the blue (~ 4260-7020A), and the VPH-Red grism (with the OG590 
filter) to cover Ha and [Nii] in the red (~ 5800-9800A). The spectral resolution was 3.1 (6.3) A 
FWHM in the blue (red). For most of our LDSS3 targets we only obtained spectra in the blue due 
to observing time constraints. For objects with double stellar components identified either from 
our NIR imaging or from SDSS images, we oriented the slit to go through the two stellar compo- 
nents; for the others we oriented the slit along the major axis of the galaxy. For a few targets whose 
[Oiii] emission along the primary slit position was either not spatially resolved or was particularly 
complex, we obtained spectra at a second slit position usually perpendicular to the primary slit 
position. We list slit positions, total exposure times, and seeing measured from field stars in ac- 
quisition images in Table [B We took wavelength calibration spectra and flat fields after observing 
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each object. During the course of each night, we observed two white dwarfs at different airmasses 
for spectrophotometric calibration and several K and M giants for velocity calibration. 

Our DIS observations were carried out during eight nights between 2009 July 15 and 2010 
July 14 UT. The observing conditions were partly cloudy or clear but not photometric on the nights 
of 2009 July 15, September 23, December 12 and 19, and 2010 July 07 and 14, and photometric on 
the nights of 2010 February 15 and June 10 UT. During most nights the seeing was poor, ranging 
between 1".5 and 2". 7 (except the nights of 2009 December 12 and 19 with seeing around 1".0). 
DIS has a 4'x6' FOV and 0".414 pixels. The spectral resolution was 1.8 (1.3) A FWHM in the 
blue (red) channel. We adopted a \".5x6' slit and the B1200+R 1200 gratings centered at 5000 and 
7000 A (or 5200 and 7050 or 5500 and 7450 or 5550 and 7550 A, depending on target redshifts). 
The slit orientation was determined in the same way as in our LDSS3 runs. 

We reduced the LDSS3 and DIS data following standard IRAlQ procedures jTodyll 19861) and 



with the COSMOS reduction pipeline@. The 2d data reduction included bias subtraction, flat field- 
ing, cosmic ray removal, wavelength calibration and spatial rectification, flux calibration and ex- 
tinction correction, sky subtraction, and aligning and stacking of individual frames. We applied 
telluric correction over the extracted Id spectra using standard stars; the O2 A-band (7580-7740A) 
is the major feature that affects the region of interest in the red. The quality of the LDSS3 data is 
substantially better than that of the DIS data due to the much better seeing and larger telescope, 
hence we will always present LDSS3 data whenever available. 



3. Interpretation 

Our combined NIR imaging and slit spectroscopy data reveal that double-peaked [Om] emis- 
sion arises from a diverse set of circumstances. We classify them by examining the NIR images 
and the two-dimensional spectra simultaneously. Given the typical seeing conditions in our NIR 
imaging, two stellar components with luminosities that are within an order of magnitude of each 
other can be identified at separations > 0.6". The seeing is poorer in our optical slit spectroscopy. 
However, the two [O iii] components are spectrally resolved, which makes it easier to deblend them 
spatially in the 2d spectrum. We measure the spatial offset between the peak emission of the two 



IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 

^http://www.ociw.edu/Code/cosmos 
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Table 1 . Observing Summary 



NIR Imaging Slit Spectroscopy 





Instrument 


Seeing 


UDS UclLc 


Exp Time 


Instrument 


Seeing 


UDS L^tlie 


PA 


Exp Time 






(, ) 


T TT 


(min) 




( ) 


T TT 


(, ) 


(sec) 


0002+0045 


PANIC K, 


0.48 


091230 


9 


DIS 


1.1 


091212 


71 


1800 


0009-0036 


PANIC K, 


0.63 


091231 


18 


LDSS3 


0.94 


100113 


278 


900 


0116-1025 


PANIC Ks 


0.52 


091229 


18 


LDSS3 


0.83 


100112 


95 


1800 


0135-0058 


PANIC Ks 


0.58 


091229 


18 


LDSS3 


0.86 


100112 


101 


1800 


0135+1435 


PANIC / 


0.55 


100101 


18 


LDSS3 


0.92 


100113 


200 


900 


0156-0007 


PANIC Ks 


0.64 


091229 


18 


LDSS3 


1.0 


100113 


299 


1800 


0400-0652 


PANIC Ks 


0.48 


091229 


18 


LDSS3 


0.58/0.57 


100112/100114 


95/173 


2700/1800 


0837+1500 


PANIC Ks 


0.66 


091229 


27 


LDSS3 


1.1 


100113 


322 


1800 


0851+1327 


PANIC Ks 


0.49 


091230 


18 


LDSS3 


0.91 


100114 


201 


900 


0942+1254 


PANIC Ks 


0.50 


091229/100102 


36 


LDSS3 


0.88/0.95 


100112/100113 


213/348 


2700/1800 


0958-0051 


PANIC Ks 


0.45 


091230 


18 


LDSS3 


0.87 


100113 


200 


900 


1009+0133 


PANIC Ks 


0.44 


091230 


18 


LDSS3 


0.84 


100114 


217 


900 


1019+0134 


PANIC Ks 


0.70 


091229 


18 


LDSS3 


0.94 


100114 


303 


900 


1038+0255 


PANIC Ks 


0.45 


091231 


18 


LDSS3 


0.75 


100114 


270 


900 


1108+0659 


PANIC Ks 


0.54 


091229 


18 


LDSS3 


0.86/0.95 


100112/100113 


320 


4800 


1131-0204 


PANIC Ks 


0.47 


091231/100102 


36 


LDSS3 


0.94/0.71 


100112/100113 


267 


2700/3600 


1146-0226 


PANIC Ks 


0.64 


100102 


18 


LDSS3 


0.98 


100114 


233'!' 


900 


1146+5110 


2MASS K 








DIS 


1.7 


100215 


51 


4800 


1322+2631 


PANIC Ks 


0.73 


100530 


18 


DIS 


2.0/1.8 


100707/100714 


79 


1800/3600 


1332+0606 


PANIC / 


0.52 


100101 


18 


LDSS3 


1.2/0.86 


100112/100114 


196 


2400/2100 


1341+2219 


PANIC K, 


0.74 


100530 


18 


DIS 


1.4 


100707 


175 


3600 


1356+1026 


PANIC a;. 


0.79 


100530 


18 


LDSS3* 










1450+0838 


PANIC Ks 


0.85 


100530 


18 


DIS 


2.4/2.5 


100707/100714 


55 


3600/1800 


1552+0433 


PANIC /f, 


0.71 


100530 


18 


DIS 


1.3 


100610 


303 


3600 


1556+0948 


PANIC Ks 


0.58 


100530 


18 


DIS 


1.5 


100610 


158 


3000 


1630+1649 


PANIC a;. 


0.86 


100530 


18 


DIS 


1.6 


100610 


39 


3000 


2252+0029 


PANIC /f, 


0.83 


100530 


27 


DIS 


1.5 


091219 


162 


2400 


2255-0812 


PANIC Ks 


0.86 


100530 


18 


DIS 


1.5 


091219 


180 


2400 


2304-0933 


PANIC a;. 


0.80 


100530 


18 


DIS 


2.0 


090923 


61 


1800 


2310-0900 


PANIC Ks 


0.62 


091230 


18 


DIS 


2.0 


090923 


73 


2700 


2333+0049 


PANIC ^, 


0.76 


091229 


18 


DIS 


1.5 


091219 


117 


3600 



Note. — Summary of our NIR imaging and optical slit spectroscopy observations. The full SDSS designations are given in Tabl e 2. For some 
object s we have two observations, which are listed separately, t Parallactic angle. *The spectroscopic observation was reported in lOreene et al J 

Hont). 
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Table 2. Object Properties 



SDSS Name 


redshitt 


PANIC olfset " (kpc) 


Spec offset " (kpc) 


Category 


000249.07+004504.8 


0.0868 




0.8 (1.3) 


NLR kinematics 


000911.58-003654.7 


0.0733 




< 0.2 (< 0.3) 


ambiguous 


011659.59-102539.1 


0.1503 




0.9 (2.3) 


NLR kinematics 


013546.93-005858.5 


0.1595 




0.2 (0.6) 


NLR kinematics 


013555.82+143529.7 


0.0719 




0.8 (1.1) 


NLR kinematics 


015605.14-000721.7 


0.0806 




0.6 (0.9) 


NLR kinematics 


040001.59-065254.1 


0.1707 




1.3/<0.2 (3.8/< 0.6) 


NLR kinematics 


083713.49+150037.2 


0.1408 




1.3 (3.2) 


NLR kinematics 


085121.94+132702.2 


0.0931 




0.6(1.0) 


NLR kinematics 


094205.83+125433.7 


0.1543 




0.2/< 0.2 (0.5/< 0.5) 


ambiguous 


095833.20-005118.6 


0.0860 




0.8 (1.3) 


NLR kinematics 


100921.26+013334.6 


0.1437 




0.2 (0.5) 


ambiguous 


101927.56+013422.5 


0.0730 




0.3 (0.4) 


ambiguous 


103850.13+025555.1 


0.0762 




0.9(1.3) 


NLR kinematics 


110851.04+065901.4 


0.1816 


0.5 (1.5) 


0.9 (2.7) 


binary AON" 


113126.08-020459.2 


0.1463 


0.6(1.5) 


0.6 (1.5) 


binary AGN" 


114610.04-022619.2 


0.1225 




1.3 (2.9) 


NLR kinematics 


114642.47+511029.6 


0.1300 


2.7 (6.2) 


2.5 (5.8) 


binary AGN'' 


132231.86+263159.1 


0.1441 




2.1 (5.3) 


ambiguous 


133226.34+060627.4 


0.2070 


1.5 (5.1) 


1.5 (5.1) 


binary AGN" 


134114.87+221957.8 


0.1152 




1.7 (3.5) 


NLR kinematics 


135646.11+102609.1 


0.1231 


1.3 (2.9) 


1.3 (2.9)* 


binary AGN* 


145050.60+083832.6 


0.1168 




< 0.4 « 0.8) 


ambiguous 


155205.93+043317.5 


0.0803 




1.2(1.8) 


NLR kinematics 


155619.30+094855.6 


0.0678 




0.4 (0.5) 


ambiguous 


163056.75+164957.2 


0.0341 




0.8 (0.5) 


NLR kinematics 


225252.94+002928.4 


0.1525 




<0.4« 1.1) 


ambiguous 


225510.12-081234.4 


0.1494 




0.4(1.0) 


ambiguous 


230442.82-093345.3 


0.0319 




0.8 (0.5) 


NLR kinematics 


231051.95-090011.9 


0.0944 




< 0.4 « 0.7) 


ambiguous 


233313.17+004911.8 


0.1699 




<0.4« 1.2) 


ambiguous 



Note. — Properties of the 31 objects in our sample. In those cases which showed two nuclei in NIR imaging, column 3 shows the spatial offset 
(in units of " and kpc) between the two nuclei, which is measured from the emission peaks of the two nuclei. Column 4 shows the spatial offset 
(in units of " and kpc) between the two velocity components of the naiTow line emission, measured from the emission peaks of the two velocity 
components in the slit spectrum. The smallest spatial offset of the two velocity components we can measure is the pixel scale, i.e., ~ 0.2" for 
LDSS3 and ~ 0.4" for DIS. For 10400-0652 and J 0942+1254 we h ave measurements for two slit position angles. The last column shows our 
classification of these objects (see "Published in lLiu et al] i2010al) ; *Slit spectrum reported in lGreene et al.l |201 ll) . 
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velocity components in the 2d spectra (Table [21), which can be measured down to one pixel scaled 
(0.188" for LDSS3 and 0.414" for DIS). 

Based on the NIR imaging and slit spectroscopy data, we classify the 31 objects in three 
categories (Table [2l): 

1. Kpc-scale binary AGNs: There are two spatially resolved (or marginally resolved) stellar nu- 
clei in NIR imaging. The two velocity components of the narrow line emission are spatially 
coincident with the pair of nuclei seen in the NIR. These are the best candidates for kpc-scale 
binary AG Ns in our sample. Five objects are included in this category, four of which were 



reported in lLiu et al.l (|2010ah . 



2. NLR kinematics in single AGNs. The NIR imaging shows a single smooth stellar compo- 
nent, but the two velocity components of the narrow lines are spatially offset by > 0.6". 
We performed tests with simulated images of two stellar components with different lumi- 
nosity contrasts, structural parameters and seeing conditions, and concluded that two stellar 
components with comparable (0. 1 < L1/L2 < 10) luminosities under the actual seeing would 
have been resolved in NIR imaging at the separation indicated by the two [Oiii] components. 
There could still be two stellar components hidden in the system if, e.g., the luminosity con- 
trast is larger than an order of magnitude, the separation of the two stellar components is 
smaller than that inferred from the narrow line emission, or one of the stellar components 
has an unusual surface brightness profile. But given the similar appearances of these objects 
in NIR imaging and slit spectroscopy to the classic example of Mrk 78 (see ®, the simplest 
explanation is that the double peaks emerge from multiple kinematic components in a single 
NLR for the majority of these sources. Fifteen objects are included in this category. We refer 
to these objects below as having NLR kinematics origin. 

3. Ambiguous cases. For the remaining objects, a single nucleus is seen in the NIR imaging, 
and the two velocity components of the narrow lines are spatially offset by < 0.4", smaller 
than the resolution of our NIR imaging. These objects could be single AGNs with NLR 
structure on smaller scales, or binary AGNs at smaller separations. We need better spatial- 
resolution observations and/or different slit positions to test these scenarios. Eleven objects 
are included in this category. 

These three categories represent our best effort to interpret these objects based on current data; 
they are by no means exact. In particular, the binary AGN and NLR kinematics classifications 



^Measuring sub-pixel offsets by fitting model line profiles to the data requires both high signal-to-noise ratio and 
symmetric spatial line profiles. These criteria are generally not satisfied by our spectroscopic data. 
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both have caveats that may cause us to misinterpret their nature. Keeping this in mind, we adopt 
this classification scheme in our following discussion. Below we summarize the objects in each 
category, and discuss the caveats. 



3.1. Kpc-Scale Binary AGNs 

Our primary interest is to identify bona fide kpc-scale binary AGNs, which motivated our 
original search for these double-peaked objects in SDSS spectroscopic database. Most of our NIR 
imaging preceded the slit spectroscopy, and we obtained spectra of all our targets with resolved 
double nuclei in the NIR. Five out of ~ 60, or ~ 10% of our objects show resolved double stellar 
nuclei in their NIR images. Their two-dimensional (2d) spectra show that the two [Qui] velocity 
components seen in the spatially-integrated spectra are spatially coincident with the two stellar 
continuum peaks. 



For the sky-subtracted NIR imaging data, we use galfit (|Peng et al.ll2010|) to model the light 
distribution with multiple components. PSFs were taken from stars within the same image. Due to 
the complexity of merging systems, and the fact that the results depend on the quality of the data, 
we only use the exponential disk (expdisk), de Vaucouleurs (devauc) and Sersic (sersic) profiles in 
each fits, and we urge caution on the interpretation of the best-fit model for some of these objects. 
During the fits we did not fix any of the model parameters, and for each fit we tried different 
combinations of the above three profiles until it reaches the minimum reduced x^. Fig.[T]shows the 
best-fit model (and the residuals) along with data for the four objects with R^NIC data. The best-fit 
models are summarized in Tabled where the best- fit parameters have been rounded using the la 
statistical uncertainties from galfit. Due to the nature of nonlinear multi-component models and 
possible systematic effects involved in the data (such as bad PSF or sky subtraction), the statistical 
uncertainties reported by galfit are an approximation of the actual uncertainties at best. We use 
these model fits to estimate the luminosity ratios of multiple stellar components in these merging 
systems, and to determine the galaxy type of each component. 

For the slit spectra, we extract one-dimensional (Id) spectra at different spatial locations 
along the slit. We model the Id spectra in different spatial bins with double-Gaussian (or double- 
Lorentzian if a better fit can be achieved) profiles for the double-peaked lines, plus a power-law 
model for the local continuum. Fig. [3] shows an example of the results of our modeling. Our goal is 
to see if there are observable trends in the velocity, line width, and line flux ratio of the two narrow 
line components. Unfortunately, given the typical seeing in our slit spectroscopy, these trends are 
only apparent in a few cases. There are several objects for which the spectral quality is too poor to 
perform such analysis. Below we briefly comment on individual objects. 
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J1108+0659. This object was reported in iLiu et al.1 (llOlOal) . It shows two nuclei in its NIR 
image, separated by ~ 0.5". The two [Oiii] velocity components are spatially coincident with the 
two NIR nuclei. Its ^^-band image and 2d spectrum for the [Oiii]-H/3 region are shown in Fig. 
I 2l The pair of n uclei are well resolved in a recent NIR adaptive-optics (AO) imaging observation 
(|Fu et al.ll2010|) . The best-fit AT^-band model consists of two de Vaucouleurs bulges embedded in 
an exponential disk. The luminosity contrast of the two bulges is ~ 0.15 magnitude (0.06 dex). 



J1131-0204. This object was reported in iLiu et al.l (|2010a|) . It shows two resolved stellar 
nuclei separated by ~ 0.6". The slit spectroscopy confirmed the coincidence of the two [Oiii] 
velocity components with the two NIR nuclei. Its ^^-band image and 2d spectrum for the [Oiii]- 
H/3 region are shown in Fig. |4l The two stellar nuclei are embedded within a galactic disk, and the 
long "spur" features seen in the 2d spectrum to > 15 kpc are ionized gas emission from the disk 
which traces the galactic rotation curve. The best-fit ^^-band model consists of one Sersic bulge 
(n ^ 3) and one de Vaucouleurs bulge embedded in an exponential disk. The luminosity contrast 
of the two bulges is ~ 2 magnitude (0.8 dex). The model fit is imperfect and there is some residual 
spiral structure indicative of interactions. 



J1146+5110. This object was reported in iLiu et al.l (|2010a|) . The double stellar nuclei were 
marginally resolved in the 2MASS image, and our slit spectroscopy subsequently confirmed the 
coincidence of the two [Oiii] velocity components with the two NIR nuclei. This is also one of 
the few cases that show resolved optical double nuclei in the SDSS images. The southern nucleus 
itself seems to have a complex NLR geometry and t wo NLR veloci ty components. The double 
nucleus was well resolved in the NIR AO imaging in IFu et al.l (|2010l) . with a luminosity contract 
of 0.7 magnitude (0.28 dex). We do not have a R^.NIC image for this object. 



J1332+0606. This object was reported in iLiu et al.l (|2010ah . The two stellar nuclei were 
clearly resolved in optical (SDSS) and 7-band images, and are spatially coincident with the two 
[Oiii] velocity components seen in our slit spectroscopy. Its 7-band image and 2d spectrum for the 
[Oiii]-H/3 region are shown in Fig. Ul The best-fit 7-band model consists of three components: a 
Sersic bulge {n ^ 2) for the southern nucleus, an exponential disk for the northern nucleus, and 
an exponential disk with a scale-length of ~ 1.6". The luminosity ratio of the three components 
is ~ 2 : 1 : 1. The less luminous northern stellar component corresponds to the stronger [Oiii] 
emission component seen in the slit spectrum, as in J 1322-1-2631. 



J1356+1026. This object was studied in lGreene et al.l (|201l|) . It shows two continuum sources 
in the optical separated by 3 kpc, correspond ing to the north and south knots of [Om] emission 
seen in the slit spectrum in iGreene et al.l (1201 ih . The tw o knots of rO i n] ern ission have a relative 
velocity offset of ~ 200 kms"^ The slit spectrum in iGreene et al.l (|201l|) shows a rich [Om] 
emission structure, including a giant [Om] bubble to the South of the southern continuum. Its 
^s-band image is shown in Fig. |9l This is a galaxy with a highly disturbed morphology and its 
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IRAS fluxes indicate that it is a ULIRG. The best-fit ^^-band model consists of three components: 
a Sersic bulge (n ^ 5) for the southern nucleus, a Sersic bulge (n ^ 3) for the northern nucleus, 
and an exponential disk component (r^ ~ 0.8") towards the north-western corner. The luminosity 
ratio of the three components is ~ 15 : 5 : 1. However, due to the highly disturbed morphology of 
this syste m, we urge cauti on on the best-fit model. This object was also observed with NIR AO 



imaging in 



Fuet all ( 20101) . 



All of these objects appear to be major mergers (the ^^-band luminosity ratio between the two 
main stellar components is less than 10), which is partly due to our selection based on comparable 
[Oiii] luminosities. But we also note that the more massive stellar component does not necessarily 
have stronger [O iii] emission. 

We classify these objects as kpc-scale binary AGNs based on spatially coincident double stel- 
lar nuclei and [Oiii] emission. While the binary scenario seems to be the most natural explanation, 
it is possible that only one SMBH is active and ionizing the gas clouds in both nuclei. In this 
single- AGN scenario, the [Oiii] gas clouds in the non-AGN host are further away from the ioniz- 
ing source than those in the AGN host. This difference will lead to different ionization states in 
the [Oiii] clouds in the two hosts. If we assume that the electron density is similar in both hosts, 
we expect ve ry different [Oiii l/H/? flux ratios of the two narrow line components, contrary to what 



we observe (LiuetaL 2010a!) . However, the electron density may well be quite different in the 



two hosts, allowing a single AGN within one host to be responsible for the [Oiii] emission in both 
hosts. This is particularly relevant for the ULIRG J1356-I-1026, where the interstellar medium 
(ISM) may be quite clumpy. Our current data are insufficient to completely rule out the single 
AGN possibility. We are currently acquiring images with HST and Chandra, as well as higher S/N 
slit spectroscopy for these objects, and further investigations of these objects will be presented in 
future work. 



3.2. NLR Kinematics in Single AGNs 

The remaining (~ 90%) objects with NIR imaging data do not show resolved (or marginally 
resolved) double nuclei at the limit of our resolution (~ 0.6"). These unresolved cases fall into one 
of the following categories: a) they are binaries at smaller projected separations; b) they are two 
accreting BHs each with its own NLR, co-rotating within a single merged stellar bulge, or c) they 
are single AGNs with complex NLR kinematics. 

About 60% of the objects that appear single in the NIR imaging show spatially resolved [Oiii] 
emission (typically > 0.6") in the 2d spectra, which correspond to the two velocity components 
seen in the spatially-integrated spectra. Double stellar nuclei separated on these scales would have 
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Table 3 . Model fits of NIR images 



Object 








q 


PA 


Mag 






(pixel/kpc) 


(pixel) 




(°) 


(Vega) 


J 1108+0659 (Ks) 






x'=0.19 








devauc (S) 


4 


3.2/1.2 


(94.1,99.8) 


0.60 


46 


15.46 


devauc (N) 


4 


3.2/1.2 


(97.3,103.6) 


0.61 


-42 


15.32 


expdisk 


1 


8.5/3.2 


(101.4,99.8) 


0.80 


46 


15.70 


T1 1 ^1 —0904 (1^ ^ 

J 1 1 J 1 U^UH- ^,A^J ) 






Xi> — 








sersic (E) 


3.2 


25.2/8.1 


(98.1,101.6) 


0.46 


-88 


15.12 


devauc (W) 


4 


14.8/4.7 


(104.2,100.3) 


0.28 


-69 


17.24 


expdisk 


1 


22.2/7.1 


(101.5,97.9) 


0.71 


13 


15.67 


J 1332+0606 (7) 






X?.= 1.08 








expdisk (N) 


1 


3.0/1.3 


(98.2,107.6) 


0.53 


-33 


17.92 


sersic (S) 


2.0 


A All. 9 


(102.2,96.4) 


0.85 


27 


17.06 


expdisk 


1 


13.1/5.5 


(100.1,115.9) 


0.55 


46 


17.79 


J1356+1026* (Ks) 






x'. = 0.91 








sersic (N) 


3.2 


4.6/1.3 


(99.0,108.0) 


0.79 


-80 


15.55 


sersic (S) 


5.6 


63.1/17.4 


(100.2,97.4) 


0.53 


-8 


14.34 


expdisk 


1 


6.7/1.9 


(115.6,121.2) 


0.67 


-1 


17.32 



Note. — Best-fit surface brightness models for four objects that show multiple stellar com- 
ponents in the NIR. Only the Sersic (including de Vaucouleurs) and the exponential disk pro- 
files were used in these fits. We report best-fit parameters for the Sersic index («), the effective 
radius (Re) or scale-length (r^), the centroid of each component (xc,yc)> the aspect ratio {q), 
the position angle of the major axis (PA), and the integrated magnitude for each component 
(normalized using the total 2MASS flux). These parameters are rounded using the statistical 
errors from the fits. Due to the complexity of these systems, we caution that a "successful" 
model may not be the unique model, and the actual errors of these parameters are expected 
to be substantially larger. This is especially a concern for J 1356 (marked with a "*"), whose 
morphology is quite disturbed. We report scales in units of pixels, where 1 pixel corresponds 
to 0.125". 
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been identified in our NIR imaging. Thus the lack of spatially coincident stellar nuclei favors either 
scenario (b) or scenario (c). However, the NLR dynamics are presumably affected by the bulge 
potential more than by the BH, thus it is difficult to maintain two distinct NLRs in scenario (b). On 
the other hand, NLR kinematics involving rotation or outflows on sub-kpc to kpc scales are quite 
common for local Seyferts (see ® . In fact, in some of the objects with good spatial quality we can 
see velocity/velocity dispersion gradients along the slit direction in the 2d spectrum (see below), 
which strongly supports the kinematics scenario. Therefore we believe that the double velocity 
peaks in the vast majority of these objects arise from complex kinematics in a single NLR . We 
now comment on each of these objects in detail. 

J0002+0045. For this object the two velocity components of [Oiii] are spatially offset by 
~ 0.8", but they do not have a corresponding pair of nuclei seen in the NIR image. Fig. [TOl shows 
the Ks image and the 2d spectrum. Note that this object has a north-east companion ~ 5" away, 
which does not have [O iii] emission; this companion galaxy is at the same redshift as J0002+0045 
measured from stellar absorption features in the slit spectrum. 

J0116-1025. This object has two spatially offset [Oiii] emission peaks, which correspond to 
the blue- and red-shifted components in the double-peaked line profile respectively. The spatial 
offset of the two [Om] components is ~ 1.1", and there is no corresponding double nucleus seen 
in the NIR image. Fig. [TT] shows the image and the 2d spectrum. The host galaxy clearly has 
a disk component. There is also a small galaxy about ~ 4" away from the center of JOl 16-1025 
to the West, which was not covered by our slit observation. Fig. [12] shows the Id slices of the 2d 
spectrum at different distances from the peak of the continuum emission. The [O iii] A5007/H/3 flux 
ratio is almost independent of position, but the [Om] line width increases towards the center of the 
continuum emission. No obvious velocity gradient is seen for either of the two [Om] components. 
Combining the NIR imaging and slit spectroscopic data, this object is best explained by a rotational 
[Om] disk co-planar with the stellar disk. 

J0135-0058. Fig.[T3lshows the image and the 2d spectrum for this object. The long slit was 
placed to cover the tidal feature to the south-east of the galaxy seen in the NIR image, and hence 
was misaligned with the major axis of the disk. Fig. [14] shows the Id slices of the 2d spectrum 
at different distances from the peak of the continuum emission. The [Om] A5007/H/3 flux ratio is 
almost constant until the outermost apertures (6 and 7), where the [Om] A5007/H/3 decreases as it 
is now tracing the faint tidal feature seen from the 2d spectrum shown in Fig. \T3\ There is a slight 
velocity gradient for the blue- and red-shifted components. Although the spatial offset between 
the two [Om] components is only ~ 0.2" (possibly due to the misaligned slit position angle), we 
believe a rotational [Om] disk is the best explanation given the velocity gradient seen in the 2d 
spectrum as well as the apparent disk morphology seen in the NIR. 

J0135+1435. Fig.[T5]shows the image and the 2d spectrum for this object. The two velocity 
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components of [Oiii] are spatially offset by ~ 0.8", while no corresponding pair of nuclei was seen 
in the NIR. Fig.[T6]shows the Id slices of the 2d spectrum at different distances from the peak of the 
continuum emission. The [Oiii] A5007/H/3 flux ratio increases towards the center of the continuum 
emission. Velocity gradients for both [Oiii] components are clearly seen in Fig. [161 and there is 
some indication of increasing line width toward the center of the continuum emission, although 
decomposition into two components is not always successful at each aperture. This object is best 
explained by a rotational [Oiii] disk, with an asymptotic flat rotation velocity VcS'mi ~ 200 kms"'. 

J0156-0007. Fig. [17] shows the image and the 2d spectrum for this object. A disk mor- 
phology is apparent in the NIR image. The two velocity components of [Oiii] are spatially offset 
by ~ 0.7", with no corresponding pair of nuclei seen in the NIR. Fig. [T8] shows the Id slices of the 
2d spectrum at different distances from the peak of the continuum emission. The [Oiii] A5007/H/3 
flux ratio is almost constant at each aperture location. Weak velocity gradients can been seen for 
both [Oiii] components. At most locations the line width is narrow, which suggests that coherent 
rotation is the dominant motion. This object is best explained by a rotational [Oiii] disk, with an 
asymptotic flat rotation velocity sin/ ~ 150kms"^. 

J0400-0652. This object shows a smooth single profile in the Kg image, while its slit spectra 
sh ow spatia l ly res olved [Om] emission extending to ~ 3". A subsequent NIR AO image obtained 



by lFu et al.l (120 10|) did not reveal a double nucleus at ~ 0. 1" resolution. In Fig. [19] we show its Kg 
image and 2d spectra at two position angles, and in Fig. [20] we show the Id spectral diagnostics. 
The 2d spectrum along R^= 95° shows rather complicated [Om] emission region kinematics. It 
shows a high velocity dispersion near the center, and some velocity gradient along the slit direc- 
tion, which is indicative of disk rotation and/or outflows. Clearly there are more than two [Om] 
components. On the other hand, the 2d spectrum along FA= 173° shows much less structure, which 
is presumably caused by the different spatial coverage. This object is best explained by NLR kine- 
matics, possibly involving both outflows and rotation. A more detailed follow-up of this object 
(i.e., with IFU or multiple-slit spectroscopy) is highly desirable to resolve the [Om] kinematics 
map. 

J0837+1500. Fig.[2T|shows the image and the 2d spectrum for this object. The two velocity 
components of [Om] are spatially offset by ~ 1.5", with no corresponding pair of nuclei seen in 
the NIR. Fig. [22] shows the Id slices of the 2d spectrum at different distances from the peak of the 
continuum emission. 



J0851+1327. Fig. [23] shows the image and the 2d spectrum for this object. The two 
velocity components of [Om] are spatially offset by ~ 0.6", with no corresponding pair of nuclei 
seen in the NIR. This galaxy has an edge-on disk component seen in the NIR, which is also seen 
in its SDSS optical image and the 2d spectrum shown in Fig. [24] This object is best explained as a 
rotational [Om] disk, given the disk morphology seen in the NIR and in the optical. 
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J0958-0051. Fig. [25] shows the image and the 2d spectrum for this object. A disk com- 
ponent is seen in the NIR image. The two velocity components of [Oiii] are spatially offset by 
~ 0.8", with no corresponding pair of nuclei seen in the NIR. Fig. [26] shows the Id slices of the 2d 
spectrum at different distances from the peak of the continuum emission. Velocity gradients are 
clearly seen for both [Oiii] components. This object is best explained by a rotational [Oiii] disk. 

J1038+0255. Fig. [27] shows the image and the 2d spectrum for this object. The two velocity 
components of [Oiii] are spatially offset by ~ 1", with no corresponding pair of nuclei seen in the 
NIR. This object is best explained by either a rotational [Oiii] disk or biconical outflows. 

J1146-0226. Fig. [29] shows the Ks image and the 2d spectrum for this object. The image 
shows a smooth single stellar bulge, while the 2d spectrum shows a clear rotation-curve like veloc- 
ity gradient across the slit. Under poor spatial resolution or if the object were observed at higher 
redshifts, the rotation-curve would be unresolvable, and the [O iii] emission would appear as two 
distinct velocity components spatially offset by ~ a few kpc. This object is best explained by a 
rotational [O iii] disk. 

J1341+2219. Fig. [3T] shows the image and the 2d spectrum for this object. A disk compo- 
nent is seen in the NIR image. The [Oiii]-H/3 part of our 2d spectrum has low quality so we show 
the Ha region instead. The two NEL velocity components are spatially offset by ~ 1.7", with no 
corresponding pair of nuclei seen in the NIR. This object is best explained by a rotational [O iii] 
disk. 

J1552+0433. For this object the two [Oiii] components are spatially offset by ~ 1.2" in 
their peak emission, as seen from its 2d spectrum shown in Fig. [32] Its image shows a disk 
component, which seems to be warped at the north-west edge. Nevertheless, two stellar bulge 
components separated by ~ 1.2" would have been seen in the NIR image if this object were a 
kpc-scale binary AGN. Thus we classify this object in the NLR kinematics category. 



J1630+1649. Fig. [33] shows the image and the 2d spectrum for this object. The two [Om] 
velocity components are spatially offset by ~ 0.8", with no corresponding pair of nuclei seen in 
the NI R. A recent NIR AO image of this object at ~ 0. 1" resolution did not show a double nucleus 



either (|Fu et al.|[2010h . This object is best explained by a rotational [Om] disk. 



J2304-0933. For this object the NIR image shows a clear disk component. The seeing was 
quite poor for the slit spectroscopy observation, but the two NEL components seem to be spatially 
offset by ~ 0.8" in their pea k emission. A recent NIR AO image of this object reveals no double 



nucleus at ~ 0. 1" resolution (|Fu et al.|[2010|) . It is thus most likely a case of NLR kinematics. 



While many of these NLR kinematics cases show a clear disk morphology in imaging data, 
a few objects, such as J0837-I-1500 and Jl 146-0226, show an early-type morphology. The ionized 
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gas emission in these early-type hosts is extended on kpc scales. Emission line kinematics studies 
of local early-type galaxies often show ionized gas em ission to kpc scal es, and the gas kinematics 
is usually decoupled from the stellar kinematics (e.g.. ISarzi et al.ll2006() . which is consistent with 
our findings here. However, the emission line strength is much stronger in our AGNs, which may 
indicate a larger gas reservoir in these early-type hosts than their local counterparts. 



3.3. Ambiguous Cases 

About 40% of objects that appear single in NIR imaging were classified as ambiguous cases, 
i.e., the combined NIR imaging and slit spectroscopic data are insufficient to distinguish between 
the binary scenario and the kinematics scenario. In essentially all of these cases, the two velocity 
components of the [Oiii] emission have undetectable or small spatial offset that is below the resolu- 
tion limit of the NIR imaging (i.e., < 0.4")- In the binary scenario, the two BHs and their associated 
stellar components would have smaller projecte d separations, w hich would need higher-resolution 



observations such as ground AO imaging (e.g.. lFu et al.ll2010l) or space-based observations to re 



solve the pair of nuclei. Nevertheless the real separation between the two BHs is unlikely to be 
much less than ~ kpc because NLRs have intrinsic sizes of hundreds of pc or above. In the case 
of NLR kinematics (rotation or outflows), either the relevant scale is below one or two kpc, or the 
slit position happens to be perpendicular to the outflow axis or the rotational disk (cf., see Fig. [19] 
for J0400-0652). Fig. |36] shows the NIR images and 2d spectra for these objects. We now briefly 
comment on individual objects in this category. 

J0009-0036. In this object the two [Oiii] components are spatially offset by < 0.2" in the 
2d spectrum, and the NIR image shows a single nucleus. Its NIR image shows a minor structure 
~ 1" south-east from the center of the peak NIR emission, which is not contributing to the [Oiii] 
emission. A better slit position and higher resolution imaging and spectroscopy are needed to rule 
out the existence of a companion at < 0.2" separation. 

J0942+1254. This object was observed with two slit positions. Both of the 2d spectra show 
no spatial offset (< 0.2") between the two [Oiii] velocity components. No double nucleus was 
seen in the NIR image. However, the optical and NIR images of this object show spectacular tidal 
features, indicative of recent mergers. This object deserves a higher-resolution observation, and 
we suspect it is likely to be a small projected separation binary AGN. 

J1009+0133. This object shows a minor companion ~ 1" north-west of the peak emission in 
the NIR image. Unfortunately due to our slit position, we couldn't confirm that these two stellar 
nuclei correspond to the two [O iii] components seen in the 2d spectrum. Although the NIR image 
alone is suggestive of a binary, a new slit position is needed to test this hypothesis. 
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J1019+0134. The NIR image shows a disturbed morphology in the outskirts, indicative of 
a recent merger event. The two [Om] components are spatially offset by ~ 0.3" as shown in the 
2d spectrum. Seeing in the NIR is ~ 0.7". This could either be a binary or NLR kinematics. 
New slit spectroscopy with a different position angle and/or better spatial resolution are needed to 
distinguish the two scenarios. 

J1322+2631. For this object we only have a slit spectrum covering the Ha region. The two 
narrow line components are spatially offset by ~ 2.1" as shown in the 2d spectrum. In the 
band image we see two main stellar nuclei separated by ~ 2.3", and there is also a third (minor) 
companion to the southwest, about 2.5" from the central brightest stellar component. A close 
examination of the 2d spectrum suggests that the two [O m] components are probably not spatially 
coincident with the two main stellar nuclei, even though their separation happens to be consistent 
with that of the two stellar nuclei. If this were true, it would imply a kinematic origin for the 
double-peaked [Oni] emission in this object. However, the seeing in our 2d optical spectrum is 
poor and the locations of the stellar continua are not well determined. Hence it is still possible that 
the two spatially offset [Oni] components are coincident with two of the three stellar components 
seen in the NIR image. The poor seeing in the spectrum also prohibited measurements of the 
absorption redshifts of the northeast and southwest companions, whose continua are blended with 
the continuum of the brighter central source. A new slit spectrum with better spatial resolution is 
needed to draw firm conclusions on this object. 

J1450+0838. The 2d spectra shows that the two [Om] components are consistent with no 
spatial offset (< 0.2"). The NIR image shows a single nucleus at ~ 0.8" resolution. 

J1556+0948. The two [Om] components are offset by ~ 0.4" in the 2d spectrum, and the 

NIR image shows a single nucleus at ~ 0.6" resolution. The double-peaked emission lines in this 
source are likely to be due to NLR kinematics, but higher-resolution NIR imaging data is needed 
to rule out the binary scenario. 

J2252+0029. The two NEL components are consistent with no spatial offset in the 2d spec- 
trum. Higher-resolution imaging observation and/or different slit positions are required to distin- 
guish the binary and the NLR kinematics scenarios. 

J2255-0812. The NIR image shows a disturbed disk component, and the two NEL compo- 
nents show a small spatial offset, < 0.5", in the 2d spectrum. Extended line emission is seen on 
opposite sides of the central continuum, which is caused by ionized gas in the rotating galactic 
disk. This object is likely a case of NLR kinematics, but a higher resolution imaging observation 
and/or different slit positions are needed to rule out the binary scenario. 

J2310-0900. This object has a companion ~ 4" away, which does not contribute to the [Om] 
emission seen in the SDSS fiber spectrum. The companion is too faint to measure its redshift in our 
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slit spectrum. The seeing was quite poor for the slit spectroscopy observation, and the two NEL 
components are consistent with no spatial offset. A better quality slit spectrum and/or different slit 
positions are needed to distinguish the binary and the NLR kinematics scenarios. 

J2333+0049. This object has a minor companion about 6" to the southeast, which is too faint 
to measure its redshift in our slit spectrum. The two NEL components are consistent with no spatial 
offset in the 2d spectrum. It appears single in our NIR imaging and in the recent NIR AO imaging 
( Fu et aLlEoiO ) at ~ 0.1" resolution. 



4. Discussion 

4.1. The Bulk Properties of Kpc-scale Binary Type 2 AGNs 

Bearing in mind the small-number statistics of our sample, in Fig. |37] we compare the bulk 
properties of kpc-scale binary AGNs (circles) and those that originate from NLR kinematics (crosses), 
measured from spatially-integrated SDSS fiber spectra (Paper I). The left panel shows the FWHM 
of each [O iii] component as a function of the fiber-integrated stellar velocity dispersion. The red 
and blue symbols are for the redshifted component and the blueshifted component, respectively; 
the gray symbols are for all the objects in our parent sample (Paper I). The five binary AGNs all 
reside in the right half of the plot, with larger stellar velocity dispersions on average. This is consis- 
tent with the existence of multiple stellar systems in these binary AGNs. The middle panel shows 
the velocity offset of the two [Om] components as a function of the stellar velocity dispersion. The 
right panel shows the total [Om] A5007 luminosity versus the [Om] velocity offset. The binary 
AGNs seem to occupy a different luminosity regime from those classified as NLR kinematics. Fig. 
|37] suggests that the bulk properties of kpc-binary AGNs might be statistically different from those 
with an origin in NLR kinematics, which, if true, can be used to refine the selection of binary 
candidates from spatially-integrated spectroscopy. However, our current sample of kpc-scale bi- 
nary AGNs is still small, and a larger statistical sample of binary AGNs is needed to confirm these 
trends. 

In these binary systems, the smaller stellar component does not necessarily correspond to the 
more luminous [Om] component. If BH mass is proportional to bulge mass/luminosity and [Om] 
luminosity is proportional to the intrinsic AGN luminosity, then the smaller BH must be accreting 
at a substantially higher Eddington ratio than the larger BH in Jl 131-0204 and J1332-I-0606. It is 
likely that the smaller component in these merging systems is more gas-rich and/or easier for the 
nuclear region to be disturbed during the merger, and hence the BH accretion is more efficient. 
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4.2. The Frequency of Kpc-scale Binary Type 2 AGNs 

Although we have only followed up a small fraction of the double peaked [O iii] type 2 AGNs 
in our sample (Paper I), the ~ 30 objects with both NIR imaging and slit spectroscopic data pre- 
sented in this paper allow us to estimate the fraction of kpc-scale AGNs among these double- 
peaked [Oiii] AGNs. For most of our targets we carried out the imaging first, and promising 
targets seen in the imaging data were preferentially observed spectroscopically. We have obtained 
slit spectra for only about half of the imaged objects. Therefore we must estimate the fraction of 
binary AGNs and NLR kinematics taking into account the spectroscopic completeness. Of the 61 
objects that we have observed in the NIR, five were classified as kpc-scale binary AGNs with re- 
solved double nuclei that are within the 3" diameter of SDSS fibers ^ and subsequently confirmed 
with slit spectroscopy. Hence the frequency of kpc-scale binary AGNs among these double-peaked 
[Oiii] objects is ~ 10%. Of the remaining ~ 90% of the objects, ~ 50% are best explained by NLR 
kinematics ( §13.21) and ~ 40% are ambiguous cases. At least some of the latter case are of NLR 
kinematics origin, and the remainder are probably binary AGNs with somewh at smaller sep ara- 



tions, and could be id entified with be tter spatial resolution observations (e.g., IFu et al.ll2010l) . If 



we take the results of iFu et al.l (120 10|) and assume all the objects with resolved double nuclei in 
their NIR AO imaging are binary AGNs (but see S4.3I) . then the binary AGN fraction among the 
double-peaked [Oiii] objects increases to ~ 20%. It is conceivable (although unlikely) that kpc- 
scale binary AGNs with even smaller separations would not have been revealed with NIR AO 
imaging, and if we assume this entire ambiguous sample is binaries, we get a hard upper limit of 
the kpc-scale binary AGN fraction of ~ 50%. Therefore, a conservative estimate of the genuine 
kpc-binary AGN fraction of the double-peaked [Oiii] type 2 AGNs is 10% < /i < 50% and the 
upper bound should be considered as a hard limit. Our observations thus indicate that the majority 
of these double-peaked [Oiii] objects reflect NLR kinematics and are not kpc-scale binary AGNs. 
However, in some of the systems that we classified as kinematic cases, we see disturbed morphol- 
ogy or close companions, which may indicate that the (single) AGN activity was triggered by a 
recent merger event. 

Only ~ 1% of low-redshift (z < 0.3) type 2 AGNs have double-peaked [Oiii] profiles in their 
spatially-integrated spectra (Paper I). If the orbit inclination is random for binary AGNs, then some 
of them will not show a double-peaked [O iii] profile since the line itself has a finite width (due to 
the intrinsic line width and instrument broadening), and thus will not enter our sample. The real 
situation is more complicated. The selection completeness depends on the spectral resolution of 
SDSS, and the intrinsic distributions of binary properties, such as the relative velocity of the two 



^Multiple stellar components separated by more than 3" is unlikely to be covered within the SDSS fiber, hence do 
not correspond to the two [O iii] velocity components seen in SDSS spectra. 
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BHs at ~ kpc separations, and the FWHM and strength of individual narrow Une components. The 
SDSS spectral resolution will exclude double-peaked [Om] objects with a relative velocity offset 
< ISOkms"', and our sample in Paper I is strongly biased against weak [Om] objects. In Paper 
I we performed simple Monte Carlo simulations of mock spectra to estimate the selection com- 
pleteness of kpc binaries with double-peaked [Om] with assumed distributions of kpc-scale binary 
AGN properties, combined with the SDSS spectral resolution and random line-of-sight (LOS). We 
found that the selection completeness rapidly decreases with decreasing the LOS velocity offset, 
decreasing the [Oin] equivalent width, and increasing the line FWHM; and all these effects are 
well coupled. 

Here we perform similar Monte Carlo simulations, with reasonable properties that are appro- 
priate for our sample. We assume a physical relative velocity offset of 500 km s"^ between the two 
velocity components (i.e., both NLRs are on circular orbits with circular velocities of 250km s"^), 
a fixed FWHM= 250 km s"^ for both components, and assign an [Om] A5007 rest equivalent width 
EW[oiii]5007 = 3EW[oiii]4959 Uniformly distributed within 15 -45 A for each component. The as- 
sumed circular velocity of the binary is reasonable given that our objects are massive galaxies with 
stellar mass ~ 10'^ Mq, and the adopted FWHM and EW values are also typical of our sample 
objects (Paper I). We further assume random orientations of the binary orbits, and generate mock 
spectra with the assumed [Om] properties, a continuum S/N~ 5pixer\ and the SDSS spectral 
resolution. We visually inspect these mock spectra and flag double-peaked objects. We found 
~ 20% of the simulated objects are identified as double-peaked [Om] objects. As discussed ear- 
lier, 10-50% of double-peaked objects have compelling evidence for two active nuclei (/i). Now 
assuming that we can detect only 20% of true double-peaked objects, we expect that 5% of the 
z < 0.3 type 2 AGN population is double, and thus that a fraction 0.5% < /2 = 5% x /i < 2.5% of 
the total population is composed of kpc-scale binaries. It is very likely that there are populations 
of kpc-scale binary AGNs with a relative (unreduced) orbital velocity < 150kms"' or with more 
discrepant [Om] luminosity ratios that are missed in our SDSS sample. 

At face value, this low fraction of kpc-scale binary AGNs is at odds with the hypothesis that 
AGNs are triggered by major mergers. It either implies that the fraction of time that a binary AGN 
spends during this kpc phase is very small compared to the total AGN lifetime, or only a tiny 
fraction of kpc-scale binary AGNs have two distinct NLRs with comparable luminosity and are lit 
up simultaneously. Alternatively, it could be that major mergers are not the dominant mechanism 
for triggering low-luminosity AGN activity at low redshift, since the fueling rate is not as stringent 
for AGNs as for luminous quasars, and alternative fueling routes (i.e., secular processes and minor 
mergers) may suffice. Refining the binary AGN fraction at various separations and combining 
with numerical simulation results are needed to sort out these possibilities, and our current study 
represents the first step towards quantifying the demographics of binary AGNs on kpc scales. 
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4.3. The Different Behaviors of Gas and Stars 



Our NIR imaging and optical slit spectroscopy demonstrate the importance of combining 
imaging and spectroscopy data in confirming the binary AGN nature of double-peaked objects. 
Objects that show spatially offset [Om] emission in slit spectra are not necessarily kpc-scale bi- 
nary AGNs, as seen in most of the cases ascribed to NLR kinematics discussed in ^3.2[ Likewise, 
objects that show apparent double nuclei in the NIR imaging are not necessarily kpc-scale binary 
AGNs, as the companion may not have [O iii] emission at all. Objects such as J0002-I-0045 pre- 
sented in this paper clearly demo nstrate the latte r case. For the same reason, not all of the objects 
with resolved multiple nuclei in IFu et al.l (|2010f) are binary AGNs, even if these multiple nuclei 
are covered in a single S PSS fiber. For example, J1157-f- 0816 shows three components in both 
optical and NIR imaging ([Smith et al.ll2010l : IFu et al.ll2010|) . However, neither the upper or lower 
companions in the slit spectrum (Fig.[38l) show [Om] emission (although both are at the redshift of 
the central object), and the double-peaked [Om] emission is likely caused by the NLR kinematics 
of the central source. 



5. Conclusions 

We have presented NIR imaging and optical slit spectroscopic data for 31 double-peaked 
[Om] type 2 AGNs drawn from the sample of Paper I. The combination of imaging and spec- 
troscopy allows us to identify the origin of the double-peaked [Om] emission in these peculiar 
objects. Our main conclusions are: 

• ~ 10% of these objects show spatially coincident double stellar nuclei and NLR emission. 
The two spatially offset [Om] emission components are also separated in velocity, and lead 
to the double-peaked profile seen in the spatially-integrated spectra. These objects are best 
explained by a pair of merging AGNs separated on kpc scales. 

• ~ 50% of these objects show a single smooth stellar component at ~ 0.6" resolution, while 
the two [Om] velocity components are spatially offset by > 0.6". For objects with sufficient 
spectral quality we usually see velocity (or velocity dispersion) gradients indicative of a 
kinematic signature. These objects are best explained by NLR kinematics in single AGNs, 
involving rotation and/or outflows. A recent study of a sin gle double-peaked [O m] AGN 



(SDSS J15 17-1-3353) also favors the kinematics explanation (|Rosario et al.ll2010l) . 



• The remaining ~ 40% of these objects show a small (< 0.4") spatial offset between the two 
[Om] velocity components and a single smoo th stellar com ponent. Some of these objects 



are likely binaries at smaller separations (e.g. jFu et al.ll2010t) . with the remainder likely due 
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to NLR kinematics at smaller scales. Follow-up observation with better spatial resolution is 
needed to distinguish between the two scenarios. 

• Our observations demonstrate the necessity of combining imaging and slit spectroscopy data 
to identify bona fide kpc-scale binary AGNs. In particular, objects that show spatially offset 
[Oiii] components are not necessarily binary AGNs, nor are those objects that show spatially 
resolved stellar components. 

• We estimate that ~ 0.5-2.5% of the local (z < 0.3) type 2 AGNs are kpc-scale binary AGNs 
in a major merger, where both components are obscured type 2 AGN with comparable [Om] 
luminosities. 

Our follow-up observations of the double-peaked [Om] sample demonstrate the feasibility 
of selecting kpc-scale binary AGN candidates using this particular spectroscopic feature, with a 
> 10% success rate. So far we have only followed up about one third of our parent sample and 
already doubled the number of known kpc-scale binary AGNs. Our imaging and spectroscopy 
follow-up is still ongoing, accompanied by detailed investigations on individual systems in other 
wavelengths and/or with better spatial resolution. A natural step forward would be AO imaging 
and IFU spectroscopy of merging systems and systems with outflows to obtain detailed spatial and 
kinematic maps. By the end of our survey, we expect to increase the number of known kpc-scale 
binary AGNs by an order of magnitude. The increased statistics of the kpc-scale binary AGN 
sample will improve the estimate of the binary fraction and the study of bulk properties of kpc- 
scale binary AGNs, and will shed light on the effects of mergers on the host and AGN activity 
in these systems. In the meantime, there is strong need to explore the parameter space in merger 
simulations to confront the observed statistics (e.g., Blecha et al. 2010, in preparation). 
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A. Model Fits for the NIR Images of Kinematic and Ambiguous Cases 

For completeness we here provide the galfit results for objects that we classified as NLR 
kinematics around single AGNs or ambiguous cases. Even though these AGNs do not appear 
to have substantial multiple components within the host galaxy, reproducing the exact surface 
brightness profile is still challenging for some cases. We restrict ourselves to a basic model of a de 
Vaucouleurs bulge + an exponential disk for the whole galaxy. We also tried a Sersic bulge plus an 
exponential disk to improve the fits if possible. In a few objects there are one or two companions 
overlapping with the light from the galaxy, and we fit additional components to these companions. 
The fitting results are shown in Figs.[39](for the NLR kinematics cases) and|40](for the ambiguous 
cases). In the NLR kinematics cases, no substantial residuals are seen at the locations indicated by 
the offset [Om] emission seen in the 2d spectra. 



This preprint was prepared with the AAS LMj5f macros v5.2. 
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Fig. 1. — Model fits of the NIR surface brightness for four objects classified as binary AGNs 
( ^3.11) using GALFiT. North is up and East is left. The first column shows the data, the second 
column shows the models and the last column shows the residuals. See the text in ^3.1l for details 
regarding the model fits. Note that the residual maps have been re-stretched to enhance the contrast 
and the residuals are not important compared to the observed fluxes (the apparent excess in some 
of the residual maps has values less than 5% of the observed fluxes at the corresponding locations). 
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Fig. 2.— SDSS Jl 108+0659 (binary). Top: PANIC NIR image in K,. North is up and East is left. 
The magenta line shows the direction of the slit used in the spectroscopy, with the arrow indicating 
the upper direction in the extracted 2d spectrum. This object has two (marginally) resolved stellar 
nuclei separated by ~ 0.5" in the cental galaxy. There is also a faint companion to the right which 
is not covered by the slit spectrum. Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with 
corresponding lines marked (note that the locations of these line marks are approximate). The two 
velocity components are spatially offset by 0.9". The stellar continua are not separated in the 2d 
spectrum due to the limited seeing and the proximity of the two stellar nuclei. 
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Fig. 3. — Diagnosis of the 2d spectrum of Jl 108+0659. The left panel shows the Id spectra summed over the entire 
slit (uppermost) and for slit slices at 7 different spatial locations (flux scale is arbitrary). The dotted vertical lines divide 
the [Oiii] A5007 line into the blueshifted and redshifted components. The dashed vertical lines separate the emission 
line regions from the continuum regions. The upper-right panel shows the spatial profiles of the continuum (black) and 
the [O III] A5007 emission line region (cyan), while the blue and red lines show the spatial profile of the blueshifted 
and redshifted [O iii] A5007 emission. We slice the 2d spectra in seven spatial bins (marked by the numbers), centered 
on the spatial peak of the continuum emission, and the Id spectra for these spatial bins are shown in the left panel. 
For the Id spectrum in each spatial bin, we try to deblend the two velocity components in the emission lines by fitting 
the spectrum with double Gaussian/Lorentzian functions and a local continuum model. We consider the deblending 
successful if the two components have a flux ratio greater than 0.3 and both components were detected at > 3a. If 
the deblending is successful we over-plot the two velocity components in blue and red respectively in the aperture 
spectra, otherwise we over-plot the whole model in cyan. We plot the velocity offset (relative to the systemic redshift) 
and FWHM of the two components, and the [O iii] X5QQ1/Hf3 flux ratio measured from the model in the bottom-right 
panels. In a few cases the apparent "successful" deblending in a spatial bin can still be spurious as judged by eye; 
we nevertheless plot results for these bins, but such bins should be ignored in assessing the spatial gradients of these 
quantities. 
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SDSSJ1131-0204 z=0.146 
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Fig. 4.— SDSS Jl 131-0204 (binary), f/pper.- PANIC NIR image in isT,. This object has two 
stellar nuclei separated by ~ 0.6". Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with 
corresponding lines marked (note that the locations of these line marks are approximate). The two 
velocity components are spatially offset by ~ 0.6". Notation is the same as Fig.[2l For this object, 
the two stellar nuclei are embedded in a disk, and the long "spur" features seen in the 2d spectrum 
are from ionized gas emission in the disk. 
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Fig. 5. — Diagnosis of the 2d spectrum of Jl 13 1-0204. Notation is the same as Fig.[3l The spatial 
distribution of the continuum shown in the upper-right panel has extended wings, which arise from 
the prominent disk as seen in the optical and NIR images. There are no coherent trends seen in the 
velocity offsets, line widths and [Oiii]/H/3 flux ratio, as functions of location. 
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Fig. 6. — SDSS Jl 146+51 10 (binary). Upper: DIS 2d spectrum for the H/3-[0iii] region. Bottom: 
DIS 2d spectrum for the Ha region with corresponding lines marked (note that the locations of 
these line marks are approximate). We do not have a PANIC image for this object, but its 2MASS 
image shows two nuclei separated by ~ 2.7". The two continua of this system are apparent in the 
bottom spectrum due to the relatively large separation of the two nuclei. The line emission in the 
southern nucleus itself has two velocity components. The line emission in the northern nucleus 
is spatially offset by ~ 2.5" from the continuum of the southern nucleus, and is very weak in the 
H/9-[0iii] region due to poor spectral quality. 
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Fig. 7.— SDSS J1332+0606 (binary). Upper: PANIC NIR image in J. This object has two 
stellar nuclei separated by ~ 1.5". Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with 
corresponding lines marked (note that the locations of these line marks are approximate). The 
two velocity components are spatially offset by ~ 1.5". Notation is the same as Fig.[2l The two 
continua are also separated in the 2d spectrum. The northern nucleus with the weaker continuum 
has stronger [Oiii] emission. 
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Fig. 8. — Diagnosis of the 2d spectrum of J1332+0606. Notation is the same as Fig. [3l The two 
velocity components are relatively well separated in space, and the stronger redshifted component 
has a larger line width than the weaker blueshifted component. 
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Fig. 9.— SDSS J1356+1026 (binary). Upper: PANIC NIR image in K^. This object has two stellar 
nuclei separated by ~ 1.3". It shows a highly disturbed r norphology, p r esuma bly caused by the 
merger. The slit spectroscopy of this object was reported in loreene et al. ( 2011 ). 
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Fig. 10.— SDSS J0002+0045 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
shows a single nucleus, and has a companion ~ 5" away toward the north-east (which is not covered 
by the SDSS fiber). Bottom: DIS 2d spectrum for the H/3-[0iii] region with corresponding lines 
marked (note that the locations of these line marks are approximate). The two velocity components 
are spatially offset by ~ 0.8". The north-east companion does not have detectable [Om] emission. 
Notation is the same as Fig. [21 



Fig. 11.— SDSS JOl 16-1025 (NLR kinematics). [/;?;?er.- PANIC NIR image in i^,. This object 
has no resolved double nucleus, and it has a small companion ~ 4" away towards the west (which 
is not covered by the SDSS fiber). Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with 
corresponding lines marked (note that the locations of these line marks are approximate). The two 
velocity components are spatially offset by ~ 0.9". Notation is the same as Fig.[2l 
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Fig. 12. — Diagnosis of the 2d spectrum of JOl 16-1025. Notation is the same as Fig. [3l The 
[Oiii]/H/3 flux ratio is almost independent of position (within ~ 0.15 dex), but the [Oiii] line 
width seems to increase towards the center of the continuum emission except for the blueshifted 
component in bins 5 and 6. No obvious velocity gradient is seen for either of the two [Oiii] 
components, and the data are consistent with a flat rotation curve with Vc sin? ~ 200 kms"^ on each 
side of the continuum. 
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Fig. 13.— SDSS J0135-0058 (NLR kinematics). Upper: PANIC NIR image in K,. This object has 
no resolved double nucleus. There are tidal features towards the south-east, which is covered by our 
slit observation. Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with corresponding lines 
marked (note that the locations of these line marks are approximate). The two velocity components 
are spatially offset by ~ 0.2". Notation is the same as Fig.[2l 
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Fig. 14. — Diagnosis of the 2d spectrum of J0135-0058. Notation is the same as Fig. [3l The 
[Oiii]/H/? flux ratio is almost constant until the outermost apertures (6 and 7), where it drops 
rapidly as it is now tracing the faint tidal feature seen in the 2d spectrum. There is a slight velocity 
gradient for the blue- and red-shifted components. Deblending of the emission lines is ambiguous 
in bins 3 and 4. 
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Fig. 15.— SDSS J0135+1435 (NLR kinematics). Upper: PANIC NIR image in K^. This ob- 
ject has no resolved double nucleus. Bottom: LDSS3 2d spectrum for the H/?-[Oiii] region with 
corresponding lines marked (note that the locations of these line marks are approximate). The 
two velocity components are spatially offset by ~ 0.8". Velocity gradients can be seen in the 2d 
spectrum. Notation is the same as Fig. |2l 
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Fig. 16. — Diagnosis of the 2d spectrum of J0135+1435. Notation is the same as Fig. [3l We 
detect clear velocity gradients for the blueshifted and redshifted velocity components, and possible 
increasing velocity dispersion towards the center of the continuum. The [Oiii]/H/3 flux ratio also 
increases towards the center of the continuum in general. 
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Fig. 17.— SDSS J0156-0007 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has a single nucleus, and it has a large disk component. Bottom: LDSS3 2d spectrum for the 
H/3-[0iii] region with corresponding lines marked (note that the locations of these line marks are 
approximate). The two velocity components are spatially offset by ~ 0.6". Notation is the same as 
Fig.H 
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Fig. 18. — Diagnosis of the 2d spectrum of J0156-0007. Notation is the same as Fig.[3l Deblend- 
ing of the two velocity components is ambiguous in several spatial bins. The [Oiii]/H/3 flux ratio 
varies within ~ 0.2 dex in bins for which we have measurements. Weak velocity gradients can 
been seen for both [Oiii] components. 
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Fig. 19.— SDSS J0400-0652 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has no resolved double nucleus. Two slit positions are shown. Bottom: LDSS3 2d spectrum for 
the H/3-[0iii] region with corresponding lines marked (note that the locations of these line marks 
are approximate). The 2d spectrum shows a rich kinematic structure at one slit position, but not at 
the other slit position. Notation is the same as Fig. |2l 
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Fig. 20. — Diagnosis of the 2d spectra of J0400-0652. Notation is the same as Fig.[3l The upper 
panel is for the slit position PA= 95° and the bottom panel is for the slit position PA= 173°. The 
two slit positions show very different distributions of the line emission. Deblending of the two 
velocity components is ambiguous in most spatial bins because of the complicated kinematics as 
seen in the 2d spectrum. 
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Fig. 20. — Continued. 
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Fig. 21.— SDSS J0837+1500 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has no resolved double nucleus. Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with cor- 
responding lines marked (note that the locations of these line marks are approximate). The two 
velocity components are spatially offset by ~ 1.3". Notation is the same as Fig.|2l 
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Fig. 22. — Diagnosis of the 2d spectrum of J0837-I-1500. Notation is the same as Fig. [3l The 
two [Oiii] components are comparable in line strength and symmetric about the center of the 
continuum. No obvious trends are seen in the velocity offset and line width of the two [Oiii] 
components. There is some indication that the [Oiii]/H/3 flux ratio increases towards the center of 
the continuum emission, which may suggest that the ionizing source is coincident with the center 
of the continuum. 
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Fig. 23.— SDSS J0851+1327 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has no resolved double nucleus. Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with cor- 
responding lines marked (note that the locations of these line marks are approximate). The two 
velocity components are spatially offset by ~ 0.6". Notation is the same as Fig.|2l 
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Fig. 24. — Diagnosis of the 2d spectrum of J085 1-1-1327. Notation is the same as Fig. [3l This 
object has a giant disk component seen in the optical image and 2d spectrum, which is also seen 
in the extended distribution of the continuum emission in the upper-right panel. No obvious trends 
are seen for the velocity offset and line with of the two [Oiii] components. The [Oiii]/H/3 flux is 
constant within ~ 0. 1 dex. 
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Fig. 25.— SDSS J0958-005 1 (NLR kinematics). Upper: PANIC NIR image in K,. This object has 
no resolved double nucleus, and it has a large disk component. Bottom: LDSS3 2d spectrum for 
the H/3-[0iii] region with corresponding lines marked (note that the locations of these line marks 
are approximate). The two velocity components are spatially offset by ~ 0.8". There are clear 
velocity gradients in the line emission. Notation is the same as Fig. |2l 
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Fig. 26. — Diagnosis of the 2d spectrum. Notation is the same as Fig. [3l Velocity gradients can 
been seen for both the blueshifted and redshifted components, which is indicative of rotation. The 
FWHM of the two velocity components show peculiar patterns, which do not seem to be caused by 
failed deblending, and is suggestive of more complicated kinematics in this system. The [Oiii]/H/3 
flux ratio seems to rise towards the outermost bins, which may by caused by density stratification. 
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Fig. 27.— SDSS J1038+0255 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has no resolved double nucleus. Bottom: LDSS3 2d spectrum for the H/3-[0iii] region with cor- 
responding lines marked (note that the locations of these line marks are approximate). The two 
velocity components are spatially offset by ~ 0.9". Notation is the same as Fig. |2l 
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Fig. 28. — Diagnosis of the 2d spectrum of J1038+0255. Notation is the same as Fig. [31 Deblend- 
ing of the two velocity components is not successful in over half the spatial bins. 
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Fig. 29.— SDSS Jl 146-0226 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has no resolved double nucleus. Bottom: LDSS3 2d spectrum for the H/5-[0iii] region with cor- 
responding lines marked (note that the locations of these line marks are approximate). This object 
shows a clear rotation curve. The two velocity components are spatially offset by ~ 1.3" measured 
from the centroids of the blue- and red-shifted emission. Notation is the same as Fig.[2l 
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Fig. 30. — Diagnosis of the 2d spectrum of Jl 146-0226. Notation is the same as Fig.[3l Deblend- 
ing is unsuccessful in the central bins as the velocity difference between the two components is 
too small there. This can be seen clearly in the 2d spectrum. The line width and the [Oiii]/H/3 flux 
ratio seem to increase towards the center of the continuum, which is consistent with the ionizing 
source being at the center of the galaxy. 
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Fig. 31.— SDSS J1341+2219 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has no resolved double nucleus, although it has a disturbed disk component. Bottom: DIS 2d 
spectrum for the Ha region with corresponding lines marked (note that the locations of these line 
marks are approximate). The DIS spectrum for the H/3-[0iii] region has poor quality. The two 
velocity components are spatially offset by ~ 1.7" measured from the peaks of the blueshifted and 
redshifted Ha emission. Notation is the same as Fig. [21 
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Fig. 32.— SDSS J1552+0433 (NLR kinematics). Upper: PANIC NIR image in K,. This object 
has no resolved double nucleus, although it has a disturbed disk component in the optical and in 
the NIR. Bottom: DIS 2d spectrum for the H/3-[0iii] region with corresponding lines marked (note 
that the locations of these line marks are approximate). The two velocity components are spatially 
offset by ~ 1.2". Notation is the same as Fig. [21 
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Fig. 33.— SDSSJ1630+1649(NLR kinematics). Upper: PANIC NIR image in i^,. This object has 
no resolved double nucleus, and it has a faint companion ~ 5" away towards the east. Middle and 
Bottom: DIS 2d spectrum for the H/3-[0iii] and Ha regions with corresponding lines marked (note 
that the locations of these line marks are approximate). The two [Oiii] components are spatially 
offset by ~ 0.8". Notation is the same as Fig.|2l 
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Fig. 34. — Diagnosis of the 2d spectrum of J1630+1649. Notation is the same as Fig.[3l Deblend- 
ing of the two velocity components is ambiguous in most spatial bins. The [Oiii]/H/? flux ratio 
varies within ~ 0.2 dex. 
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Fig. 35.— SDSS J2304-0933 (NLR kinematics). Upper: PANIC NIR image in K,. This object has 
no resolved double nucleus, and a large disk component is apparent. Bottom: DIS 2d spectrum for 
the H/3-[0iii] region with corresponding lines marked (note that the locations of these line marks 
are approximate). The two velocity components are spatially offset by ~ 0.8". Notation is the 
same as Fig.|2l 
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Fig. 36. — PANIC NIR images and 2d spectra for objects classified as ambiguous cases. In all 
objects except for SDSS J1322+2631, the spatial offset between the two emission line velocity 
components is smaller than the typical resolution of our NIR imaging. 
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Fig. 36. — Continued. 
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Fig. 37. — Bulk properties of the five objects classified as binary AGNs (circles) compared with 
those classified as NLR kinematics in singl e AGNs (crosses) . In all panels, the gray points are for 
all objects in the double-peaked sample of iLiu et al.l (l2010b|) . Left: Distributions in the plane of 
FWHM and stellar velocity dispersion a*. The stellar velocity dispersions are for the whole galaxy 
within the 3" diameter fiber and the FWHMs are for the two velocity components (the blueshifted 
component is i ndicated in blue a nd the redshifted component is indicated in red). Measurements 
are taken from Liu et al. ( 2010b ). Middle: Distributions in the plane of a* versus velocity offset 
between the two components, Voff- Right: Distributions in the plane of the total [Oiii] A5007 
luminosity vs Voff- 
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Fig. 38. — A slit s pectrum for Jl 157+0816 , which shows three nuclei in the optical SDSS image 



(|Smith et al.ll2010h and in the NIR image (|Fu et al.ll2010h . The slit was placed across the three 
nuclei. Only the central nucleus has [Oiii] emission, which is responsible for the double-peaked 
[Oiii] profile seen in the SDSS spectrum. 
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Fig. 39. — Model fits of the NIR surface brightness for objects classified as NLR kinematics of single AGNs 03.2l i 
using GALFiT. The first column shows the data, the second column shows the models and the last column shows the 
residuals. Note that the residual maps have been re-stretched to enhance the contrast and the residuals are not important 
compared to the observed fluxes (the apparent excess seen in some of the residual maps has values less than 5% of the 
observed fluxes at the corresponding locations). 




Fig. 39. — Continued. 




Fig. 39. — Continued. 
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Fig. 40. — Model fits of the NIR surface brightness for objects classified as ambiguous cases ( ^3. 3b using galfit. 
The first column shows the data, the second column shows the models and the last column shows the residuals. Note 
that the residual maps have been re-stretched to enhance the contrast and the residuals are not important compared to 
the observed fluxes (the apparent excess seen in some of the residual maps has values less than 5% of the observed 
fluxes at the corresponding locations). 




Fig. 40. — Continued. 



